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We present the analysis of the muon events with all muon multiplicities collected during 21804 
hours of operation of the first LVD tower. The measured depth - angular distribution of muon 
intensities has been used to obtain the normalization factor, A, the power index, 7, of the primary 
all-nucleon spectrum and the ratio, R c , of prompt muon flux to that of n— mesons - the main 
parameters which determine the spectrum of cosmic ray muons at the sea level. The value of 
7 = 2.77 ± 0.05 (68% C.L.) and R c < 2.0 • 10" 3 (95% C.L.) have been obtained. The upper limit 
to the prompt muon flux favours the models of charm production based on QGSM and the dual 
parton model. 

PACS numbers: 13.85.Tp, 13.85.-t, 96.40.Tv, 96.40.-z 



I. INTRODUCTION 

The depth - angular distribution of muon intensity measured in an underground experiment is closely related to 
the muon energy spectrum at surface. Assuming the muon survival probabilities are well known for every depth and 
every muon energy at surface, the analysis of the measured depth - zenith angle distribution of intensity allows us to 
evaluate the parameters of the muon spectrum at the sea level, i.e. the normalization constant, the power index of the 
primary all-nucleon spectrum, 7, and the prompt muon flux from the decay of charmed particles produced together 
with pions and kaons in the high-energy hadron-nucleus interactions. 



Among these characteristics the value of prompt muon flux attracts a particular interest. It can be evaluated from 
the zenith-angle distributions of muon intensities, measured at various muon energies or various depths. The fraction 
of prompt muons cannot be estimated from the muon energy spectrum or depth-intensity curve measured at one 
zenith angle because the same effect can be produced cither by the prompt muons, or the decrease of 7, or both. 

The charmed particles are produced together with pions and kaons in the collisions of primary cosmic rays with 
air nuclei. They have such short live times that they decay immediately (if their energy is less than 1000 TeV) into 
muons and other particles. Thus, for them there is no competition between interaction and decay, and the prompt 
muon energy spectrum has almost the same slope as the primary spectrum. Due to the rise of the charm production 
cross-section in the energy range 100-1000 TeV, the power index of the prompt muon spectrum, j c , can be little 
lower than 7. However, possible scaling violation in the fragmentation region can increase the value of j c . Due to 
the absence of the competition between interaction and decay of charmed particles the zenith-angle distribution of 
prompt muons is almost flat, comparing with the sec9 - distribution of the conventional muons (from the decay of 
pions and kaons). This allows to estimate the fraction of prompt muons by analysing the zenith-angle distribution of 
muon intensities. 

Numerous calculations of the prompt muon flux were done (see, for example, ||l|-|8[). Different models give the 
prompt muon fluxes which vary by 2 orders of magnitude. This is due to the uncertainties in the charm production 
cross-section, er c , x-distribution of charmed particles (x = E c /Eq), produced in pA-collisions, and the branching 
ratio of charmed particle decay into muons. The most uncertain parameter, that results in the large dispersion of 
the predicted prompt muon flux, is the ^-distribution of produced charmed particles in the fragmentation region, 
important for the charm-produced cosmic-ray muons. This distribution at high energies cannot be measured precisely 
at accelerators which give the information only about small x. Thus, to check the models of the charm production, 
the experiments with cosmic-ray muons at high energies are useful. 

The search for the prompt muon flux was done with several detectors located at the surface and underground (see, 
for example, P-Q2|). In practice, it is convenient to express the prompt muon flux in terms of the ratio, R c , of prompt 
muon flux to that of pions at vertical. Since the slope of the prompt muon spectrum is close to that of pion spectrum, 
the ratio R c is almost constant for all muon energies available in the existing experiments. The experimental data, 
collected up to now, show a large variation of R c (from to 4 • 10~ 3 ). 

In a previous paper |13|] we have presented our measurement of the single muon 'depth - vertical intensity' curve 
and the evaluation of the power index of the meson spectrum in the atmosphere using the 'depth - vertical intensity' 
relation for single muons. Here we present the analysis of all muon sample which include the muon events with all 
multiplicities. The muon survival probabilities, used to obtain the value of 7 in p3| , have been presented in Q. They 
have been calculated using the muon interaction cross-sections from |lq-jl^]. After the publication of these results, 
new calculation of the cross-section of muon bremsstrahlung and of the corrections to the knock-on electron production 
cross-section have been done |18J . In the present analysis we have taken into account the corrections proposed in |lq | 
and we have estimated the uncertainties of 7 due to the uncertainties of the cross-sections used to simulate the muon 
transport through the rock. In this paper we present a more detailed evaluation of the characteristics of the muon 
spectrum at the sea level, including the ratio of the prompt muon flux to that of pions, using the depth - zenith 
angle distributions of muon intensities (I M (x, 9)) measured with LVD in the underground Gran Sasso Laboratory. The 
analysis is based on an increased statistics comparing with the previous publications. The 'depth - vertical intensity' 
relation for all muon sample and its analysis are presented in a separate paper |19| ] . 

In Section 2 the detector and the procedure of data processing are briefly described. In Section 3 the results of the 
analysis of the muon intensity distribution (I^x, 9)) are presented. In Section 4 we discuss our results in comparison 
with the data of other experiments and theoretical expectations. Section 5 contains the conclusions. 

II. LVD AND DATA PROCESSING 

The LVD (Large Volume Detector) experiment is located in the underground Gran Sasso Laboratory at a minimal 
depth of about 3000 hg/cm 2 . The LVD will consist of 5 towers. The 1st tower is running since June, 1992, and the 
2nd one - since June, 1994. The data presented here were collected with the 1st LVD tower during 21804 hours of 
live time. 

The 1st LVD tower contains 38 identical modules |2(J. Each module consists of 8 scintillation counters and 4 
layers of limited streamer tubes (tracking detector) attached to the bottom and to one vertical side of the supporting 
structure. A detailed description of the detector was given in (g0|. One LVD tower has the dimensions of 13 x 6.3 x 12 
m 3 . 

The LVD measures the atmospheric muon intensities from 3000 hg/cm 2 to more than 12000 hg/cm 2 (which corre- 
spond to the median muon energies at the sea level from 1.5 TeV to 40 TeV) at the zenith angles from 0° to 90° (on 



the average, the larger depths correspond to higher zenith angles) . 

We have used in the analysis the muon events with all multiplicities, as well as the sample of single muons. Our 
basic results have been obtained with all muon sample. This sample contains about 2 millions of reconstructed muon 
tracks. 

The acceptances for each angular bin have been calculated using the simulation of muons passing through LVD 
taking into account muon interactions with the detector materials and the detector response. The acceptances for 
both single and multiple muons were assumed to be the same. 

As a result of the data processing the angular distribution of the number of detected muons N^((j), cos 6) has been 
obtained. The angular bin width 1° x 0.01 has been used. The analysis refers to the angular bins for which the 
efficiency of the muon detection and track reconstruction is greater than 0.03. We have excluded from the analysis 
the angular bins with a large variation of depth. 

The measured A^(</>, cos #)-distribution has been converted to the depth - angular distribution of muon intensities, 
I fJ- (x,cos9), using the formula: 



I (x co6) = Hj N n(xm((l>j), cos 0j) 

^ m ' ' J2j( A ( x rn((l>j), COS e^eiXrn^j), COS 9i) • fiy ■ T) 



(1) 



where the summing up has been done over all angles <f>j contributing to the depth x m ; A(x m (4>j), cos 9i) is the 
cross-section of the detector in the plane perpendicular to the muon track at the angles (<pj, cos 9i); e(x m {<f>j), cos 9i) 
is the efficiency of muon detection and reconstruction; fiy is the solid angle for the angular bin, and T is the live 
time. We have chosen the depth bin width increasing with the depth to have comparable statistics at all depth bins 
from 3 to 10 km w.e.. Thus, the depth bin width increases from about 100 m w.e. at 3000 m w.e. to more than 500 m 
w.e. at about 10000 m w.e. The muon intensities have been converted to the middle points of the depth bins taking 
into account the predicted depth - intensity relations for different zenith angles (we have used the parameters of the 
muon spectrum at sea level which fit well the 'depth - vertical muon intensity' relation measured by LVD Q). The 
angular bin width has been taken equal to A(cos#) = 0.025. The conversion to the middle points of the angular bins 
has been done according to the predicted angular dependence for muons from pion and kaon decay. However, due to 
the small angular bins this conversion does not change angular distributions. 

III. ANALYSIS OF THE DEPTH - ZENITH ANGLE DISTRIBUTION OF MUON INTENSITY 

MEASURED BY LVD 

The data analysis has included the procedure of fitting of the measured depth - zenith angle distribution of muon 
intensity with the distributions calculated using the known muon survival probabilites (see J13| , [l4[ , and references 
therein) modified for a new muon bremsstrahlung cross-section |lq ] and muon spectrum at sea level with three free 
parameters: normalization constant, A, power index of primary all-nucleon spectrum, 7, and the ratio of prompt 
muons to pions, R c . The depth - angular distributions of muon intensity have been calculated using the equation: 

I f ,(x 1 cos9) = / P(E ll0 ,x) ■ — — -% dE^Q, (2) 

Jo ac/iO 

where P(E fl Q^x) is the probability for muon with an initial energy E^q at sea level to survive at the depth x in 
Gran Sasso rock, and — " d ^°' c ° s is the muon spectrum at sea level which has been taken according to [gl[ : 
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where the values of cos9 have been substituted by cos9* which have been taken from either |22|] or a simple 
consideration of the curvature of the Earth atmosphere. In a search for a small contribution of prompt muons it is 
necessary to know precisely the angular dependence of conventional (from pion and kaon decay) muon intensity at 
all energies of interest. In p2] cos9* — E^ r K (cos9 = l)/E° r K (cos0), where E° r K are the critical energies of pions and 
kaons. cos9* can be understood also as the cosine of zenith angle of muon direction at the height of muon production. 
The height of muon production increases from 17 km at cos9 = 1 to about 32 km at cos9 = 0. We have found that 
the values of cos9* depend on the model of the atmosphere in the range of cos9 = — 0.3. In Figure [l] we present the 



predicted angular dependences of conventional muon intensities at the energy of 10 TeV. As can be seen, all curves 
almost coincide at cosQ = 0.3 — 1. However, there is a large spread of functions at cosQ = — 0.3. The calculations 
using eq. (0) with cosO* from (22] (upper solid curve) or the treatment of the Earth curvature with a muon production 
height of 32 km (dash-dotted curve), as well as the results of [£3| (dashed curve) give quite similar results at all cos9, 
while the original calculations of p3,p3 (lower solid curve) and the treatment of the Earth curvature with a muon 
production height of 17 km (dotted curve) are far below or above other curves at small cosQ. To be independent of 
the model we have restricted the range of cosO used in the analysis to 0.3 - 1. This increases the statistical error 
of the results decreasing at the same time the systematical uncertainty related to model used. This also reduces the 
sensitivity of the experiment to small values of R c . We note that the uncertanties in the rock thickness and rock 
density are high enough at small cosQ. Moreover, large derivative of the column density with angle together with 
muon scattering effect lead to the high uncertainties of the muon flux. This also justifies our decision to restrict the 
range of zenith angles used in the analysis. 

We have added to the original formula of J2l| the term R c , which is the ratio of prompt muons to pions. Here it 
has been assumed that the power index of the prompt muon spectrum is equal to that of primary spectrum. Really, 
due to rapid rise of charm production cross-section and the possible scaling violation in the fragmentation region, the 
prompt muon spectrum may have the power index, j c , different from 7. But the value of j c depends on the model of 
charm production. To be independent of the models we have used at the first approximation the assumption: 7 C =7- 
The full formula has been multiplied by the additional normalization constant A which has been considered as a free 
parameter together with 7 and R c . 

As aresult of the fitting procedure we have obtained the values of the free parameters: A = 1.84±0.31, 7 = 2.77±0.02 
and the upper limit on R c < 2 • 10~ 3 . Here and hereafter we present the errors at 68% confidence level (C.L.) and 
the upper limits at 95% C.L. The value of \ 2 is equal to 316.7 for 330 degrees of freedom. The estimates of the 
parameters A and 7 are strongly correlated. The larger the value of 7 is, the larger the normalization factor A should 
be. Figure 3 shows the contour plot of allowed region in A — 7 - plane. The dependence of \ 2 on R c is presented in 
Figure || which was used to obtain an upper limit on R c . The errors of the parameters include both statistical and 
systematic uncertainties. The latter one takes into account the possible uncertainties in the depth and local density, 
but does not take into account the uncertainty in the cross-sections used to simulate the muon transport through the 
rock. If we add the uncertainty in the muon interaction cross-sections, the error of 7 will increase from 0.02 to 0.05 
(for the discussion about the uncertainty due to different cross-sections see p5[). This uncertainty, however, does not 
influence the upper limit on R c . We note that the energy in eq. (|3|) is expressed in GeV and the intensity is expressed 
in cm -2 s _1 sr _1 . If we restrict our analysis to the depth range 5-10 km w.e., we obtain the following values of 
parameters: A = 1.6±g;|, 7 = 2.76 ± 0.06 and R c < 3 • 10~ 3 . 

The angular distributions of muon intensities for depth ranges of interest are presented in Figure H together with 
calculations with R c = (best fit - solid curve) and R c — 2- 10~ 3 (upper limit - dashed curve). The normalizations of 
both calculations have been done independently using the fitting procedure. The data at all zenith angles are shown 
but the analysis was restricted to the range 0.3 < cosO < 1. The error bars show both statistical and systematic 
uncertainties. The calculated distributions have been obtained using the eq. (|J) and the values of cosQ* from P2fl , As 
can be seen from Figure 0, there is no evident increase of the deflection of the data points from the best fit predictions 
(R c = 0) with the increase of depth at large cos 6* as it should be if the significant prompt muon flux is present. The 
deepest depth bin is the exception. However, due to small statistics, the data at very large depth do not affect much 
the total value of \ 2 ■ 

If the formula from p3 is used for the muon spectrum at sea level instead of eq. (0), the best fit values of 7 will be 



decreased by 0.04-0.05 and will be in agreement with the previously published values for single muons 1 13 14J] analysed 
using the formula from p4[ . This difference, being comparable with our total error, is due to the factor which is 
present in the formula from |24| ] and takes into account the rise of hadron-nucleus cross-section at high energies. This 
factor appears in the calculation J24|] if the rise of the total hadron-nucleus cross-section with energy is due to the rise 
of the differential cross-section in the central region, while the scaling is conserved in the fragmentation region. This 
factor makes the muon energy spectrum steeper and the difference in the power index of muon spectrum is about 
0.04-0.05. 

Similar analysis performed for single muons shows no evidence for prompt muon flux, too. We found the same 
values of power index and upper limit to the prompt muon flux, while the absolute intensity is 10% smaller. 

IV. DISCUSSION 

From the analysis of the depth-angular and depth distributions of muon intensities measured by LVD the following 
estimates of the parameters of the muon spectrum at the sea level have been obtained: A = 1.84±0.31, 7 = 2.77±0.02 



(68% C.L.), R c < 2 • 10~ 3 (95% C.L.). The errors include both statistical and systematic errors with the systematic 
error dominating. The systematic error takes into account the possible uncertainties in the depth and local density, 
which have been estimated from the difference between the measured and predicted intensities for all angular and 
depth bins. The uncertainties of rock thickness and local density both result in the uncertainty of the column density 
and, hence, in the uncertainty of the muon flux. The distribution of fractional differences between measured and 
predicted intensities has been found to be close to gaussian with a standard deviation of about 0.04. This value 
has been assumed as a systematic error of muon intensity due to the column density uncertainty. This value is 
equivalent to the column density error of about 1% at a depth of 3 km w.e. It is obvious that the systematic error 
is more important at small depth where the statistics is high and statistical error is negligibly small. An additional 
systematic error due to the uncertainties of the cross-sections of muon interactions used to simulate the muon survival 
probabilities should be included. According to the discussion in p5f we estimate the total uncertainty in 7 as 0.05 
and in A as 0.5. The uncertainty in the cross-sections, however, does not affect the upper limit to R c . To check 
this we have fitted LVD data with the intensities calculated with muon bremsstrahlung cross-section from p5| and 
obtained the following results: A = 1.86 ± 0.32, 7 = 2.78 ± 0.02 (68% C.L.), R c < 2 • 10~ 3 (95% C.L.). The muon 
bremsstrahlung cross-section from [|15j is a little smaller than that from |lq] . This makes the muon 'depth- intensity' 
curve (with fixed A, 7 and R c ) flatter. This is compensated in the data analysis by the increase of 7. But the 
shape of the calculated angular distribution of muon intensities at any fixed depth, used to extract the value of R c , 
is not changed and, hence, the limit on ratio of prompt muon flux to that of pions remains unchanged. However, the 
absolute value of prompt muon flux (or its limit) varies with the muon cross-sections used since the flux depends also 
on normalization constant, A, and power index, 7 (see eq. (||)). 

The value of 7 obtained with LVD data is in reasonable agreement with the results of many other surface and 
underground experiments (see, for example, ]l2| . p6| pl| ) . However, the results obtained in the experiments which used 
the indirect method of the measurement of the muon spectrum, in particular, the measurement of the depth-intensity 
curve, are strongly affected by the muon interaction cross-sections and the algorithm applied to calculate the muon 
intensities. We have used the most accurate cross-sections, known at present, and the algorithm which allows us 
to calculate the muon intensities with an accuracy of 1% for a given set of muon interaction cross-sections and for 
homogeneous medium. The algorithm can influence strongly the calculated muon intensities and, then, the final 
results (for a discussion sec, for example [|25|). Thus, the observed agreement (or disagreement) in the value of 7 does 
not mean the agreement (or disagreement) in the data themselves. 

The conservative upper limit to the fraction of prompt muons, obtained with the LVD data (R c < 2 • 10~ 3 ), even 
in the simple assumption that the power index of the prompt muon spectrum, j c , is equal to that of primaries, 7, 
rules out many models of the prompt muon production, which predict a fraction of prompt muons more than 2 • 10 . 
To make this conclusion more reliable we have carried out the analysis of the depth - angular distribution of muon 
intensity using the prompt muon spectra predicted by different models (without a constant term R c ). We conclude 
that the LVD data contradict the predictions of model 1 |l| , model II [0 and model A . The predictions of the model 
3 Q, model I ||, models B, C|, recombination quark-parton model (RQPM) || and model by Q are comparable 
with the LVD upper limit, and these models cannot be ruled out. At the same time the LVD result favours the models 
of charm production based on QGSM (see, for example, |6j) and the dual parton model H, which predict low prompt 
muon flux. 

The upper limit (95% C.L.) obtained with the LVD data is lower than the value of R c found in the MSU experiment 
(R c = (2.6 ± 0.8) • 10~ 3 at £V,o = 5 TeV fl2[ ). The LVD upper limit does not contradict the values of prompt muon 
flux, obtained in Baksan flit] and KGF P] underground experiments. Our result agrees with that of NUSEX JTn] 
which did not reveal any deviation from the angular distribution expected for conventional muons. 

We point out that the LVD sensitivity to the prompt muon flux is restricted mainly by the systematic uncertainties 
connected with the uncertainties of the slant depth and local density fluctuations and the differences in the theoretical 
shape of the muon underground intensities. 

V. CONCLUSIONS 

The analysis of the depth-angular distribution of muon intensity measured by LVD in the depth range 3000-10000 
hg/cm 2 has been done. The parameters of the muon energy spectrum at the sea level have been obtained (see eq. 
(f|): A = 1.8 ± 0.5, 7 = 2.77 ± 0.05 and R c < 2 ■ 10~ 3 (95% C.L.). The errors include both statistical and systematic 
uncertainties. The upper limit to the fraction of prompt muons, R Cl favours the models of charm production based 
on QGSM H and the dual parton model [J7|, and it rules out several models which predict a high prompt muon flux. 
Similar analysis performed for single muon events revealed the same values of power index and upper limit to the 
fraction of prompt muons, while the normalization constant is 10% smaller. 
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FIG. 1. Ratio of muon intensity at cos 9 to that at cos# = 1 for 10 TeV muons at sea level versus cosine of zenith angle 
cos# calculated using different formulae: dotted curve - eq. (3) with cos#* from Earth curvature with scale height of 17 km; 
dash-dotted curve - eq. (3) with cos 9* from Earth curvature with scale height of 32 km; upper solid curve - eq. (3) with cos 9* 
from [21]; dashed curve - calculations of [22]; lower solid curve - original formula from [21,23]. 
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FIG. 2. Contour plot of allowed region in A — 7 - plane showing strong correlation between the parameters. 
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FIG. 3. Dependence of \ 2 on the ratio of prompt muon flux to that of pions. 
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FIG. 4. The dependence of the muon intensity on the zenith angle for the depth bins of most interest in the analysis and for 
all zenith angles. The data have been converted to the middle points of depth and angular bins. Solid curve - calculation with 
7 = 2.77 and R c = (best fit to the LVD data in the whole depth range, see eq. (3)); dashed curve - calculation with 7 = 2.77 
and R c = 2 • 1CP 3 (LVD upper limit). The absolute normalization of both sets of calculations has been done independently 
using the fitting procedure. The error bars include both statistical and systematic uncertainties. 



10 



